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The soft character of organic materials leads to strong coupling between molecular nuclear 
and electronic dynamics. This coupling opens the way to control charge transport in organic 
electronic devices by exciting molecular vibrational motions. However, despite encouraging 
theoretical predictions, experimental realization of such control has remained elusive. Here we 
demonstrate experimentally that photoconductivity in a model organic optoelectronic device can 
be controlled by the selective excitation of molecular vibrations. Using an ultrafast infrared laser 
source to create a coherent superposition of vibrational motions in a pentacene/C60 photoresistor, 
we observe that excitation of certain modes in the 1500 — 1700 cm -1 region leads to photocurrent 
enhancement. Excited vibrations affect predominantly trapped carriers. The effect depends on the 
nature of the vibration and its mode-specific character can be well described by the vibrational 
modulation of intermolecular electronic couplings. Vibrational control thus presents a new tool for 
studying electron-phonon coupling and charge dynamics in (bio)molecular materials. 
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The soft character of organic materials strongly influ¬ 
ences their electronic functionality. mm In these systems 
charge hopping and electronic delocalization are deter¬ 
mined by the overlap of the molecular orbitals and, there¬ 
fore, is highly sensitive to minor changes in molecular 
geometry. Hence, the electronic properties of organic ma¬ 
terials are largely determined by the interplay between the 
electronic and nuclear dynamics of the molecules, referred 
to as vibronic coupling phenomena. A growing number of 
interdisciplinary studies show that vibronic effects lie at 
the heart of a diverse class of effects in physics, chemistry 
and biology - from non-linear behavior of molecular junc¬ 
tions [21 to photophysics of vision [3] and even olfactory 
reception. [4] Vibrational motions have been postulated 
to regulate the interaction between different molecular 
electronic states by modulating inter- and intra-molecular 
couplings, by donating or accepting extra energy quanta 
HE], and by suppressing |6] or promoting [7J quantum 
interference phenomena. 

Vibronic effects were also shown to be fundamentally 
important for the conductivity of organic materials. Vi¬ 
brational motions influence intermolecular electron tun¬ 
neling probabilities m and govern a variety of non¬ 
equilibrium phenomena such as local heating nn> switch¬ 
ing [2], hysteresis, and electronic decoherence [SJQjj. This 
makes vibrational excitation a promising tool for spec¬ 


troscopy of molecular junctions HUE], tracking charge 
transfer processes in organic- and bio-electronic systems, 
and, more generally, for the development of electronic de¬ 
vices. For example, remarkable opportunities for organic 
electronics would arise from the possibility to control 
charge transport, and, thus, affect device performance by 
coherently driving nuclear motions along a pre-selected 
reaction coordinate trajectory. However, despite many 
encouraging theoretical predictions mu the experi¬ 
mental realization of vibrationally controlled electronics 
is still elusive due to the complexity of selective control 
of nuclear motions in an actual electronic junction. 

Until now, vibration-associated charge dynamics in or¬ 
ganic electronic devices have been only controlled with 
approaches that do not include mode selectivity. For 
example, the density and the equilibrium population of 
vibrational states have been varied via chemical synthesis 
of molecules with different bond structures m and via 
thermal population of low-frequency vibrations. [6] How¬ 
ever, in principle, it should be possible to access particular 
non-equilibrium nuclear or vibronic states by using instru¬ 
mentation of optical time-resolved techniques, like visible 
pump-probe 00 HZ 1 , time-resolved stimulated/impulsive 
Raman □SHE], or transient IR absorption [20]. For exam¬ 
ple, for inorganic perovskite materials [21] and molecular 
Mott insulators [22] it has been reported that selective 
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IR excitation can lead to strong modulation of the elec¬ 
tronic properties by inducing a lattice phase transition. 
Sophisticated all-optical two-dimensional (2D) photon 
echo techniques are even capable of guiding a molecu¬ 
lar system through a desired quantum superposition of 
vibronic/vibrational states [7J 1231426] . Although such 
spectroscopic methods provide a comprehensive approach 
for probing and controlling molecular motions and have 
been applied to model systems such as molecular thin 
films or solutions, they have not yet been employed to 
control functional electronic (nano)devices. 

In this work, we combine device characterization and 
ultrafast-spectroscopy methods to experimentally demon¬ 
strate that the performance of an organic optoelectronic 
system can be controlled by selectively exciting vibra¬ 
tional modes of the molecules involved in charge trans¬ 
port. As model system we use pentacene/C60 bi-layer 
photoresistors. Our experimental approach is based on 
the interferometric extension of the pump-push photocur¬ 
rent (PPP) technique. [27] In a PPP experiment, an 
optoelectronic device is illuminated by a sequence of laser 
pulses interacting with the active material in the device. 
The result of these interactions is detected by observing 
the variations in the current flow through the device as 
a function of time delay T between the pump and push 
pulses and their spectra. Thus, PPP combines the sen¬ 
sitivity and device relevance of electronic methods with 
the excitation selectivity and ultrafast time resolution of 
optical methods. Since its introduction [28, 29 , PPP has 
been applied and discussed in the context of photovoltaics 
mm, nanoelectronics, spectroscopy m , microscopy 
[32] . and molecular junction research m- In this work, 
we extend the PPP method, using the recent progress 
in ultrafast interferometry [26, [33] that allows for a pre¬ 
cise control over the time/frequency-domain structure of 
the IR optical pulses. We apply a sequence of ultrafast 
mid-IR laser pulses to create a coherent superposition 
of molecular vibrational motions inside the active layer 
of a device and correlate this excitation with the device 
performance. 

Figure[l]a)-c) describes the organic bilayer photoresistor 
model system. The active layer of the device consists 
of polycrystalline pentacene (70 nm) and fullerene C60 
(15 nm) films (Fig. [l] a)), thermally evaporated on top of 
3-, 5- or 10 pm spaced electrodes arranged in a comb-like 
geometry on a SiCU substrate (Fig. |T] b)). We chose this 
geometry rather than a sandwich-like structure, typical 
for photodiodes or solar cells, to improve the access of 
mid-IR pump-pulses to the active layer. Adding the C60 
layer was critical to enhance the photocarrier generation 
in the film. [34] 

Figure [l] d) compares the absorption spectra of pen¬ 
tacene and C60 in the IR vibrational fingerprint region 
and in the region of the optical electronic transitions. C60 
shows several distinct vibrational modes at 1180 cm -1 , 
1430 cm -1 , and 1540 cm -1 and has a comparably low op- 
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FIG. 1. The molecular electronic device characterisation: (a) 
Molecular arrangement of molecules in the pentacene crystal 
and C60 fullerene structure, (b-c) Layout and microscope 
image of the device, (d) IR absorption in the vibrational 
fingerprint region and optical absorption spectra of pentacene 
and C60. The yellow shaded contour shows a typical laser 
spectrum used for IR push, (e) Photocurrent from the device 
as a function of visible light modulation frequency; the line is 
a Cole-Cole fit with a 2.2 ms lifetime constant and a = 0.55 
dispersion parameter. 


tical density in the visible. Pentacene has a rich spectrum 
of vibrational lines in the IR and also shows strong exci- 
tonic absorption features at frequencies above 14500 cm -1 
(690nm). According to density functional theory cal¬ 
culations, the strong IR peaks at 1300 and 1345 cm -1 
are mostly associated with C=C stretching vibrations 
along the short axis of pentacene, while the weaker high- 
frequency vibrations correspond to atomic motions mostly 
aligned with the long axis of the molecule (see SI). 

The dark I-V curves of the devices are symmetric and 
roughly linear, indicating good hole injection from the 
gold electrodes to the pentacene layer (see SI). Upon 
exposure to visible light, the current flow through the 
devices strongly increases (^ 3 times under 10mW/cm 2 
illumination). Devices without a C60 layer demonstrated 
only negligible photoconductivity, which indicates that 
singlet (and triplet) [34] excitons generated after pen¬ 
tacene excitation are dissociating at the pentacene/C60 
interface and that the charge generation proceeds through 
the interfacial charge transfer states m ■ Due to the large 
electron injection barrier at the pentacene/Au interface, 
the dark current is mostly provided by holes, while under 
illumination both holes and electrons contribute to the 
photocurrent. Unlike in a typical solar cell, both elec- 
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trodes are placed below the pentacene films. Therefore, 
electrons and holes have to pass through the pentacene, 
which is known to lead to extremely long (up to seconds) 
extraction times of electrons residing in low-lying trap 
states in pentacene. [36] This notion is confirmed by the 
dependence of the photocurrent on the light-modulation 
frequency (Fig. [l] e)). Cole-Cole analysis of this depen¬ 
dence shows a typical time constant > 2 ms, which we 
interpret as the lifetime of long-lived electronic charge 
carriers. 

Figure [2] a) shows the layout of the experiment, de¬ 
signed to observe the effect of molecular vibrations on 
the charge transport through the device. The setup com¬ 
bines a 1 kHz visible-infrared ultrafast spectrometer and 
a lock-in current probe station wired to the device un¬ 
der ^ 5V external bias. First, a visible (15000cm -1 ; 
665 nm; 1.9 eV) pump pulse illuminates the device. The 
absorption of the pump light in pentacene leads to the 
build-up of excitons and charge carriers in the active layer. 
The generated carriers produce a sequence of 1 ms-spaced 
current pulses in the measurement circuit with an av¬ 
erage photocurrent J « 10 nA, detected by the lock-in 
amplifier at 1 kHz. We note that at such low current 
densities a charge-induced phase transition m can be 
excluded. The device is irradiated with a push pulse at 
certain delay times T before or after the pump pulse. The 
push pulse can promote < 1% of the molecules to the 
excited vibrational state and can also excite low-frequency 
charge-associated IR electronic transitions. [38 The effect 
of IR light on the charge separation and transport was 
detected via the variation of device photocurrent SJ. 

Figure [2] c) presents a typical PPP transient, measured 
with a single-pulse push (one interferometer arm blocked) 
at 1250 — 1500 cm -1 . When the pump was blocked we 
observed no signal due to the push only. At negative delay 
time T, when the push pulse arrives before the pump, we 
already observe a substantial increase of the current due 
to IR excitation (i.e., SJ > 0). We associate this response 
with the excitation of long-lived photocarriers that were 
generated by the preceding pump pulse that arrives ^ 
1 ms earlier. This observation is in line with the long 
collection times of trapped carriers observed for electrons 
in pentacene. 36 At delay time T — 0, the PPP response 
promptly increases as the concentration of charges in the 
cell rises due to the arrival of the new pump pulse and 
the IR push influences their dynamics. The rapid rise is 
followed by a ^ 100 ps decay component that we assign 
to the geminate recombination of newly generated charge 
pairs, which are likely to form electrostatically bound 
charge-transfer excitons. [27 

In a broadband experiment using a single-pulse push, it 
is not possible to distinguish the effects of low-frequency 
electronic excitations from the vibronic phenomena as¬ 
sociated with the interference between the molecular vi¬ 
brational motions and charge dynamics. To separate 
and address these phenomena individually, we performed 


push-frequency resolved measurements by exploiting the 
ultrafast interferometry approach. [26] Using a Mach- 
Zehnder scheme (Fig. [ 2 ] a)), the push beam is split into 
two pulses displaced in time by an interferometric delay r. 
This leads to the formation of a 1/r periodic modulation 
in the total push spectrum (Fig. 2b) which allows for 
selective excitation of different coherent superpositions of 
modes within the bandwidth of the IR light. In a typical 
experiment, for a certain pump-push delay T, the signal 
SJ/J is detected as a function of interferometric delay 
r (Fig. ! d)). The obtained interferogram is Fourier- 
transformed along the r axis to yield the action spectrum 
of the push effect. 

Figure [ 2 ] e) shows a typical frequency-resolved PPP re¬ 
sponse of a pentacene/C60 device at negative and positive 
pump-push delay times T, and with no pump (dark). At 
both delays the response consists of a number of narrow 
peaks on top of a broad featureless response, roughly 
following the IR source spectrum. We associate the broad 
feature with intraband electronic and polaronic absorp¬ 
tion, which typically spreads between 1000 and 5000 cm -1 . 
[38] The intraband excitation brings the associated charge 
carriers to a higher-lying delocalized state, thereby en¬ 
hancing their mobility, decreasing their recombination, 
and thus increasing the current output. [27] The narrow 
features in the PPP signal match well with the absorp¬ 
tion peaks of the vibrational modes of pentacene and 
C60. Therefore these features in the frequency-resolved 
PPP response are assigned to the excitation of molecular 
vibrations that modulate the electronic dynamics. In¬ 
terestingly, the broad electronic response dominates the 
PPP signal when the push arrives after the pump, while 
the vibrational features have similar amplitudes (within 
the experimental accuracy) at positive and negative T 
delays. This observation indicates that the IR electronic 
excitation substantially promotes charge separation at the 
pentacene/C60 interface, soon after exciton generation. 
At the same time, the effect of vibrational excitation is 
present for long-lived trapped charge carriers and, there¬ 
fore, does not influence charge separation, but only carrier 
de-trapping dynamics. [27^ 

We now focus on the analysis of the vibrational features 
only. The effect of broadband electronic IR excitation 
on charge dynamics in organic semiconductors has been 
investigated previously [27] . and is outside the scope of 
this paper. To study the effect of vibrational excitation 
for a broader set of vibrational modes, we use a wide push 
spectral window of 1150 — 1700 cm -1 and long r-scanning 
to obtain high frequency resolution. We also applied 
time-domain filtering (see SI) to suppress broad features 
due to electronic excitation and non-linear field-induced 
tunneling currents. 

Figure [ 3 ] a) presents the vibration-associated PPP spec¬ 
trum covering most of the IR fingerprint frequency range. 
The amplitude of the PPP response was normalized to 
the spectral density of the push pulse to allow for a direct 
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FIG. 2. The results of pump-push photocurrent measurements: (a) The layout of the time- and frequency-resolved PPP 
experiment, (b) Time- and frequency-domain representations of the IR interferometric pulse pair, matching molecular vibrations, 
(c) Broadband PPP transient for pentacene/C60 photoresistor, (d) Typical PPP interferogram. (e) Frequency-resolved PPP 
signals for pentacene/C60 photoresistor, measured at different pump-push delay times T, and a no-pump dark measurement. 


comparison of the different vibrational lines. The spec¬ 
trum is the result of several measurements with different 
push frequencies spliced together to match the amplitude 
of the effect for the 1430 cm -1 feature, which was present 
in all measurements. In the 1150 — 1700 cm -1 region, we 
observe twelve PPP peaks at frequencies that match well 
the IR-active vibrational modes of pentacene and C60. 

We note that, while the vibrations of charged pen¬ 
tacene may differ from those of the neutral molecules [39^, 
these differences should not be observed in the PPP data. 
Firstly, the shift in frequency for most individual modes 
is small 40 , [4F and, for most modes, below our frequency 
resolution 10cm -1 ). This conclusion is also supported 
by DFT calculations, see Suppl. Information, Tables SI 
and S2. Secondly, the minor shifts of the vibrational levels 
lead to a highly efficient vibrational energy transfer [42] 
between neutral and charged molecules, which allows the 
vibrational excitation of neutral pentacene to be deliv¬ 
ered to the charge trapping sites. Thirdly, according to 
Miller-Abrahams formalism, when a carrier hops from a 
radical to a neutral state all vibrational modes coupled to 
these electronic states contribute to the transfer rate. [43] 
Therefore, it is not surprising that the lattice vibrations, 
i.e. those of the neutral pentacene, are observed in the 
detrapping dynamics. 

We observe that the amplitude of the PPP response 
does not follow the intensity of the IR absorption. For 
example, the band at 1345 cm -1 possesses a much stronger 


IR absorption than the 1630 cm -1 vibration, but shows a 
weaker PPP response. This result shows that the observed 
PPP response cannot be explained by the equilibration of 
vibrational energy between modes and average heating of 
the device active layer, thus illustrating the mode-selective 
character of the PPP response. This example illustrates 
that different atomic motions couple differently to the 
charge dynamics of the system. To exclude that the non¬ 
scaling of the PPP response with IR absorption is merely 
an effect of a different orientation of the vibration dipoles 
with respect to the exciting IR light, we also performed 
angle-dependent IR absorption measurements. These 
measurements showed that the modes exhibiting very 
different PPP effect, e.g. at 1345 cm -1 and 1630 cm -1 , 
have similar dipole orientations (see SI), which rules out 
orientation effects. 

Figure |3] b) compares the effect of vibrational excitation 
on the device photoconductivity for different vibrational 
modes, obtained by normalizing the PPP response to the 
number of photons absorbed by the vibrational mode. 
In accordance with Figure 3a, the 1300 and 1345 cm -1 
modes show the weakest coupling. The higher-frequency 
vibrations of pentacene show a 5-8 times higher effect 
on the photoconductivity. For two of the fullerene vibra¬ 
tions the effect is similar to that of the high-frequency 
vibrational modes of pentacene. 

These results can be rationalized in the framework 
of the phonon-assisted Miller-Abrahams (MA) theory. 
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FIG. 3. Experimental evaluation of vibrational control effect: 
(a) The vibrational part of the PPP response, measured at 
negative delay time (—100ps). The signal amplitude is nor¬ 
malized to the spectral density of the IR push source. The 
spectrum is obtained using two PPP spectra, each covering a 
different but overlapping wavenumber range; these are scaled 
to match the amplitude of the 1430 cm -1 mode that is present 
in both spectra. For comparison, the absorption spectrum of 
the pentacene/C60 layer is presented in red. (b) The influence 
of different vibrations on device photocurrent, estimated by 
normalizing the amplitude of the PPP signal to the absorbed 
IR intensity. The change of photocurrent absolute value cor¬ 
responds to a flat 0.5J/cm 2 per cm -1 spectral density of 
exciting IR light, fully absorbed by the vibrations. The error 
bars are calculated from standard deviations for measurements 
on different devices; the number of measurements was 10 for 
1300 — 1450 cm -1 modes and 4 for all other modes. 


[43] According to this model, carrier hopping from a 
trapping state to higher-energy (more conducting) states 
takes place via absorption of a phonon with energy A 
to compensate for the energy difference between initial 
and final electronic states; the hopping rate k is defined 
by the electron-vibrational coupling constant (y) and the 
occupation number (n a) of the absorbed phonon, i.e. 

k oc v 2 n A . (1) 

In thermal equilibrium, the occupation number of a high- 
energy molecular vibration is very small: 


( 2 ) 


For a comprehensive description of the PPP response 
the interaction with the IR photons should be included 
into the MA model. However, at the conceptual level 
the effect can be understood by assuming that an IR 
excitation creates a non-equilibrium population of the 
molecular vibrational manifold; therefore, an increase in 
the hopping probability is expected. The mode-selective 
character of the PPP response is therefore defined by the 
electron-vibration coupling constants. 

To achieve mechanistic insight into the observed phe¬ 
nomena we performed a theoretical analysis of the cou¬ 
pling between the different molecular vibrations and the 
charge carriers (holes) in pentacene. In molecular systems 
these couplings can be divided into two types, i.e. local 
(Holstein-type) and non-local (Peierls-type). [44] The 
Holstein electron-phonon interaction originates from the 
modulation of the site energies by the vibrations. Only 
totally symmetric molecular vibration modes can con¬ 
tribute to this interaction. For centrosymmetric molecules 
like pentacene, the symmetric modes are not IR-active 
and will not absorb IR photons, so that we can rule out 
Holstein electron-phonon coupling effects in the PPP re¬ 
sponse. The Peierls-type electron-phonon couplings are 
associated with the dependence of the transfer integrals 
on the distances between adjacent molecules and their 
relative orientations. mi For this type of coupling, there 
are no symmetry restrictions. 

Based on previous studies [45] , we used a triclinic poly¬ 
morph [46] to represent the pentacene layer structure 
in the calculations. Figure [4] a) shows the simulated IR 
spectra for a single pentacene molecule and for the crystal 
in comparison to the experimental IR absorption. The 
agreement between the experimental and calculated vibra¬ 
tional frequencies (with typical discrepancies < 10 cm -1 ) 
allows the assignment of the vibrational modes observed 
in the frequency-resolved PPP experiment. The 1300 
and 1345 cm -1 features in the experimental spectrum are 
associated (Fig. [ 4 ] c) and SI) with in-plane ring stretching 
modes along the short axis of pentacene [47, 48j, while the 
IR peaks at 1540 cm -1 and 1630 cm -1 are associated with 
molecular deformations along the long axis of pentacene 
(Fig. [I] d) and SI). 

The non-local hole-vibration couplings are defined as 
the derivatives of the charge transfer integrals with respect 
to the vibrational coordinates, b>i = dti/dQj , and can be 
computed numerically. [ni mi Both the transfer inte¬ 
grals and electron-vibration couplings have been derived 
in a one-electron approximation (see SI for details). Our 
results indicate that there are two main transfer integrals 
contributing to charge transfer in the pentacene crystal: 
t\ = 75meV and £2 = 32meV; both are associated with 
intermolecular interactions along the herringbone direc¬ 
tions (see the red and blue arrows in Fig. |4]b). The 
other two transfer integrals oriented along the a-axis are 
substantially smaller and do not demonstrate substantial 
modulation by IR-active modes. (See SI) The derived 


n A = exp (—A /k B T). 
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FIG. 4. (a) Simulated IR spectrum for the pentacene molecule (red) and crystal (blue), superimposed on the experimental 

spectrum for comparison. (b)Molecular crystal structure and transfer integrals ti addressed in the calculations. (c,d) Eigendis- 
placements for the modes with large transition dipoles at 1288 cm -1 and 1632 cm -1 . (e) Non-local electron-phonon coupling 
constants for two largest transfer integrals, (f) Calculated vibration-induced hopping rates for the different modes that were 
addressed experimentally, compared to the experimentally observed effect of vibrational excitation on the photocurrent, SJ/ J. 


coupling constants of the IR-active modes are shown in 
Figure 4e. The couplings in the 1400 — 1650 cm -1 range 
are about 2 to 5 times larger than in the 1200 —1400 cm -1 
range. To link the variations in coupling and the probabil¬ 
ity to de-trap a charge with a vibrational excitation, we 
estimated the rates of vibration-induced charge hopping. 
In the case of two pathways, the rate is defined within 
perturbation theory as: 

ki oc (uli + i4 ti + + v'li), (3) 

where and z/G correspond to quasi-degenerate molec¬ 
ular vibrations of similar frequencies. Figure [2] f) presents 
these hopping rates together with the experimental ob¬ 
servations from Figure [3] b). The theoretical results are 
fully consistent with the experimental data. In particular, 
they capture well the mode-selective character of the phe¬ 
nomena, with the modes below 1430 cm -1 calculated to 


have a much smaller impact on charge hopping than the 
higher-frequency ones. Based on the calculations, the in- 
termolecular electronic couplings and charge transport in 
pentacene crystals are seen to be most sensitive to stretch¬ 
ing deformations along the long molecular axis, while the 
stretching deformations along the short molecular axis 
are less important. 

In conclusion, we demonstrated that the vibrational 
coupling phenomena, which play an essential role in 
molecular-scale charge transport, can be explored and 
put to action by combining optical and electronic tech¬ 
niques. Both the experiment and theoretical calcula¬ 
tions demonstrate that different non-equilibrium geome¬ 
tries and atomic motions have different effects on the 
charge dynamics. Specifically, our results show that vi¬ 
brations along the long axis of pentacene molecules lead 
to a stronger increase of hopping transport via charge 
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de-trapping than vibrations along the short axis. The 
mode-selective vibrational control of charge dynamics 
introduced here opens up a plethora of opportunities 
for basic research, including the development of high- 
mobility organic semiconductors, and the utilization of 
vibronic phenomena for ultrafast switching of organic de¬ 
vices. In addition, the mode-selective and local nature of 
our method might be particularly useful for the identifi¬ 
cation of charge transport mechanisms and pathways in 
(bio)molecular junctions. [50j 
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